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Abstract 

By  mounting  small  disc  reflectors  that  are  moveable  relative 
to  the  Inner  reflector  surface  of  a  parabolic  dish  antenna, 
nulls  can  be  generated  in  the  side  lobe  region  of  the  power 
radiation  pattern  with  minimal  distortion  effects  to  the  main 
beam.  A  physical  optics  model  of  this  antenna  system  is 
developed  to  Investigate  in  a  simplified  direct  manner  the 
phenomena  of  phase  nulling  caused  by  disc  movement. 

Array  theory  using  isotropic  radiators  is  used  to  sample  the 
aperture  distribution  to  approximate  the  far  field  electric 
field  of  the  dish.  A  physical  optics  approximation  for 
scattering  off  a  flat  metal  disc  is  used  for  discs  and  feed 
blockage  effects.  The  array  theory  term  plus  the  feed  blockage 
term  yields  the  field  intensity  quiescent  patterns,  phase  and 
amplitude,  proportional  to  |lQulescentCe,'>>)|  .  The  physical 
optics  term  for  the  discs  yields  the  field  intensity 
cancellation  patterns,  phase  and  amplitude,  proportional  to 
I^Cancellationt®  •  Th®  quiescent  patterns  are  combined 
with  the  cancellation  patterns  to  produce  the  relative  power 
pattern,  nulled  pattern,  proportional  to  llTotal^® >®)|  No 
secondary  effects  such  as  diffraction  or  edge  illumination  are 
considered. 


A  computer  code  was  written  to  implement  this  approach  and 
the  theoretical  patterns  produced.  Actual  measured  patterns  are 
compared  to  the  theoretical  patterns.  There  is  good  agreement 
between  theory  and  actual  measurements.  Finally  applications  of 
this  antenna  system  are  presented. 


This  type  of  reflector  antenna  was  originally  proposed  for  space 
applications  (Lang,  1978).  The  feedback  circuitry  of  (EMR) 
could  be  used  to  systematically  distort  the  reflector  surface  to 
form  nulls  in  selected  directions  of  the  side  lobe  region  of  the 
antenna.  A  distortion  of  the  reflector  surface  would  change  the 
phase  of  a  corresponding  region  of  the  aperture  field 
distribution  thus  generating  a  perturbation  in  the  radiation 
pattern.  By  appropriate  control  of  the  feedback  circuitry  a 
null  could  be  "moved"  in  the  desired  direction  to  eliminate  any 
interference  from  that  direction.  The  theoretical  basis  for 
this  scheme  was  investigated  by  Havens  (Havens,  1983).  It  was 
found  that  the  control  problem  of  maintaining  a  perfect  shape, 
such  as  a  paraboloid  was  extremely  difficult.  Thus,  the  use  of 
selective  distortion  appeared  too  complicated  for  practical  use. 

Historically  adaptive  nulling  in  the  side  lobe  region  has 
been  accomplished  by  arrays.  Computer  controlled,  these 
antennas  are  capable  of  steering  their  main  beams  and  nulls. 
Array  antennas,  however,  are  complex,  costly,  and  heavy  in 
comparison  with  reflector  antennas. 

"In  the  summer  of  1983  Jacavanco  demonstrated  it  was 
possible  to  produce  nulls  in  the  pattern  of  a  rigid  reflector  by 
mounting  discs  on  the  reflector  dish  and  adjusting  their 
distance  off  the  dish  by  computer  algorithm  controlled  motors  as 
shown  in  Figure  2"  (Rudlslll,  1984b: 1). 

These  antennas  are  much  simpler  to  construct  and  offer  high 
maximum  gain  per  unit  weight.  This  approach  offers  flexibility 


and  improved  signal  to  noise  ratios  in  an  already  widely  used 


Figure  2.  Computer  Controlled  Discs - Parabolic  Reflector 

(Adaptive  Reflectors,  1984c;  Research,  1985bc) 


From  January  to  November  1984,  Rudisill  developed  a 
theoretical  background  for  the  system  of  Figure  2.  He  used 
aperture  integration  and  modeled  the  disc’s  effect  on 
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the  aperture  field  as  a  phase  shift  across  a  circular  projection 
of  the  disc  into  the  aperture  plane.  The  feed  was  assumed  to 
produce  a  parabolically  tapered  electric  field.  This  approach 
required  large  amounts  of  computer  time  for  numerical 
integration  and  theoretical  results  from  the  model  compared 
poorly  to  experimental  results.  He  concluded  this  to  be  a 
result  of  too  idealistic  an  antenna  model  rather  than  error  in 
basic  theory  (Rudisill,  1984a: 7, 59). 

During  this  same  period  Jacavanco  modeled  the  system  of 
Figure  2  using  linear  array  theory.  His  model  showed  with 
reasonable  accuracy  the  salient  features  of  nulling  with  two 
discs  (Jacavanco,  1984b). 

Approach 

The  present  approach  is  similar  to  Jacavanco’ s  linear  array 
model.  It  uses  841  isotropic  radiators  spaced  0.25  wavelengths 
apart  to  simulate  the  actual  normalized  aperture  distribution, 
approximating  the  far  zone  electric  field  of  the  reflector.  The 
disc  and  feed  blockage  effects  are  approximated  using  a  physical 
optics  method  of  scattering  by  flat  metal  discs,  assumed  to  be 
perfectly  conductive  (see  Figure  3).  All  three  terms  are 
combined  to  produce  a  system  relative  power  pattern  in  the 
H-plane.  The  Relative  Power  Pattern  is  proportional  to 

IIT0T(^  »^=0)  =  IEd(®  .^=0)+Es12(®.<’=0)''EFB(®>*=0)|^  (1-1) 

(Total)  (Dish)  (Scattering  (Feed 

Discs  1&2)  Blockage) 


where:  1tot(®»^=0)  is  “the  ■total  complex  vector  electric 

field  for  the  H-plane. 


ia  the  total  complex  vector  electric 
field  from  the  reflector  without 
discs  and  feed  blockage. 


lS12t®»^=0)  is  i'*'®  complex  vector  electric  field 
from  scattering  off  the  two  discs. 


is  the  complex  vector  electric  field 
from  scattering  off  the  parasitic 
reflector. 


Figure  3.  Two  Discs  on  a  Parabolic  Reflector 

"The  effect  of  the  aperture  blocking  can  be  approximated  by 
subtracting  the  radiation  pattern  produced  by  the  obstacle  from 
the  radiation  pattern  of  the  undisturbed  apertut'e"  (Skolnick, 
1962:269).  A  description  of  the  system  follows  this  section. 
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supporting  rods  that  move  them,  relative  to  the  reflector  Inner 
surface,  lie  on  a  common  line  with  the  dish  vertex  and 
perpendicular  to  the  dipole  direction  (see  Figure  5). 


Figure  5.  Disc  Location:  Front  View  (Not  to  Scale) 

gaisxratl,Bg  s  Hailed  BslatlYS  ggysr  Pattern 

The  antenna  system  is  mounted  on  a  pedestal  that  can  move  in 
azimuth  or  elevation  planes  (see  Figure  6).  A  pyramidal  horn, 
175  feet  away  acts  as  a  transmitter/ Jammer.  The  transmitter 
illuminates  the  frontal  aspect  of  the  dish  at  3.12  6hz  and  is 
assumed  to  be  in  the  far  field  of  the  antenna.  It  is  assumed 


Figure  6.  Antenna  Mounted  on  Pedestal 
(Research,  1985bc) 

plane  waves  of  one  milliwatt  lntensitjr>  vertically  polarised  are 
intercepted  by  the  dish.  The  pedestal  is  fixed  in  the  elevation 
plane  so  that  only  azimuthal  rotation  occurs.  Since  the 
incident  E-field  is  vertically  polarized,  this  rotation  will 
generate  on  H-plane  pattern.  The  discs  are  held  flush  on  the 
inner  reflector  surface  (see  Figure  6)  and  a  system  quiescent 
power  pattern  is  recorded.  A  Scientif ic-Atlanta  Series  1580 
Pattern  Recorder  was  used.  Once  the  quiescent  pattern  is 
recorded  it  is  rewound  so  the  nulled  pattern  can  be  recorded 
over  it  in  a  different  color  ink  for  comparison  purposes.  To 
accomplish  the  nulled  pattern  the  antenna  is  rotated  so  the 
signal  is  received  at  some  arbitrary  point  in  the  side  lobe 
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region.  It  is  held  there  while  the  system  computer  finds  the 
optimum  position  for  the  discs  to  minimize  the  received  signal. 
With  the  discs  in  the  optimum  position  a  new  power  pattern 
(nulled  pattern)  is  recorded  over  the  quiescent  power  pattern. 

The  euitenna  system  uses  an  algorithm  to  minimize  the 
received  signal.  It  varies  the  discs,  one  at  a  time,  from  flush 
to  their  maximum  out  extension  (1/2  Revolution  =  1/32  inch  out 
from  reflector  surface)  through  thirty-two  possible  positions. 


i 


At  each  position  the  actual  distances  of  the  discs  from  the 
reflector  surface  are  recorded  as  electric  phase  angles  and  the 


Power  Total  (PTOT)  is  recorded  corresponding  to  the  electric 
phase  angles.  A  gradient  search  of  the  data  is  performed  for  the 
combination  of  phase  angles  and  PTOT  that  yields  the  minimum 
PTOT.  The  computer  then  positions  the  discs  via  their  motors  to 
the  optimum  position. 

The  next  chapter  will  discuss  the  theory  involved  with  the 
nulling  phenomena  caused  by  disc  movement. 
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II.  Ibspcy 


Phase  MXlXns  With  Discs 

The  tactical  situation  this  system  could  be  used  in  is  this: 


I  FRIENDLY  TRANSMITTER _ | 

Figure  7.  Tactical  Situation  for  System  Use 

The  enemy  transmitter  is  jamming  us  through  the  receiver 
antenna’s  side  lobes.  The  receiver  antenna’s  power  pattern  by 
reciprocity  is  the  same  as  its  power  pattern  when  it  transmits. 
If  a  null  can  be  shifted  in  the  receiver  antenna  power  pattern, 
in  the  jammer  transmitter  direction,  we  can  negate  the  jamming 
interference.  If  the  jamming  interference  can  be  split  into  two 
electric  field  components  of  the  same  £unplitude  but  with  phases 
180°  opposite,  they  will  cancel  each  other  completely  due  to 
the  destructive  interference  of  two  coherent  signals,  equation 
(2-1),  (two  sinusoids  at  the  same  frequency)  (Research, 1985b) . 
Refer  to  equation  (2-1)  and  Figure  8  on  the  next  page. 


Interfering  source.  The  correct  size  and  disc  location  has  been 
determined  by  Jacavanco  (Jacavanco,  1984b: 1,2).  “If  nulling  is 
desired  in  a  uniformly  illuminated  reflector  in  the  region  of 
the  pattern  where  the  response  is  down  30dB  from  the  main  beam 
peak,  then  the  total  area  of  the  reflector  minus  the  total  area 
of  the  two  discs  should  also  be  down  by  30dB.  In  an  antenna  in 
which  the  feed  horn  imposes  a  tapered  illumination  across  the 
dish  for  side  lobe  control,  the  size  and  position  of  the  discs 
should  be  further  modified.  For  the  same  performance,  a  larger 
disc  is  required  if  positioned  in  the  region  of  lower  intensity 
than  a  disc  positioned  in  a  region  of  higher  intensity  (near  the 
center)."  Another  factor  is  that  disc  size  and  placement  on  the 
reflector  surface  should  not  degrade  the  main  beam.  This  has 
been  experimentally  proven  by  Jacavanco  (Jacavanco,  1984ab).  A 
simplified  picture  of  the  entire  system  is  presented  in  Figure  9 
(only  one  disc  shown  for  clarity). 


where:  Prsq  =  Power  received 
from  signal 
thru  dish 

PrSD  ■  Power  received 
from  signal 
thru  discs 

pRjQ  =  Power  received 
from  Jammer 
thru  dish 

PrjD  =  Power  received 
from  Jammer 
thru  discs 


Figure  9.  Simplified  System 


So  the  Total  Power  received  is 


PtoT  =  PrSQ  +  PrSD  +  PrJQ  *  PRJD  (2-3) 

Therefore  the  discs  are  moved  to  minimize  PtoT  phase 
cancellation  between  Prjq  and  Prjd*  Since  Prsq  ^  Prsd 
the  signal  received  will  be  Prsq.  the  desired  signal. 

gzstem  Coordinate  Reference 

The  reference  system  used  with  the  isotropic  radiators  in 
the  aperture  plane  (xy  plane  at  z  =  q)  is  shown  on  the  next  page 
in  Figure  10.  The  z  =  0  phase  reference  is  at  the  actual  pivot 
point  for  the  antenna  system  mounted  on  the  pedestal  (see  Figure 
6).  Where  they  are  more  convenient,  spherical  coordinates  are 
used  and  coordinates  that  refer  to  the  aperture  plane  (source) 
will  be  primed  any  time  there  is  any  confusion.  Far  field 
patterns  are  presented  in  rectangular  coordinates.  Complex 
scalars  are  designated  by  underbars  and  complex  vectors  are 
designated  by  both  overbars  and  underbars.  Unit  vectors  are 
designated  by  "hats",  in  other  words,  is  the  unit  vector  in 


the  x-direction. 


Arfaz  Ibecry;  E-fieid  qi  the 


Quiz 


To  restate  the  problem:  How  does  the  phase  of  the  electric 
fields  from  the  discs  Interact  with  the  phase  of  the  electric 
fields  in  the  side  lobe  region  of  the  dish?  This  causes 
cancellation,  which  ultimately  perturbs  the  power  radiation 
pattern;  shifting  side  lobe  nulls  to  new  angles.  Reviewing  the 
approach  from  Eq  (1-1) 

IEtOT(®*<>’=0)|  ^  =  I  ED(fl><*>=0)+Esi2(9*'*»=0)-EFB(»  .<*»=0)|  ^ 


The  development  of  Ej)( 0,0=0)  follows.  The  aperture  plane 
currents  of  the  dish  alone  will  be  sampled  using  linear  array 
theory  starting  from  the  complex  vector  potential.  Ed(^»^=0) 
will  be  derived  from  the  complex  vector  potential.  The  complex 
vector  potential  is  given  by  (Stutzman  and  Thiele,  1962:455): 


// 

Source 


e ds 
4xR 


(2-4) 


where:  7  =  currents  in  aperture  plane 

R  =  distance  between  source  points  and 
observation  points 

d  =  zjl. 

\ 


Assuming  parallel  rays,  R  for  the  phase  is  determined  by 
geometrical  techniques.  R  in  the  denominator  is  approximated  by 
r,  an  average  distance  from  the  phase  reference  point  (z  =  0)  to 
the  observation  point  (Stutsman  and  Thiele,  1962:22-24). 


Therefore  the  complex  vector  potential  in  the  far  field  is 


Aff(T')  =  ±  (r) fjf  J( 


tff')  ^ 


where:  ±  (r)  =  e~J^^ 

4irr 

_/ 

r  =  vector  from  phase  reference  point 
to  source 

r  =  distance  from  phase  reference  point 
to  observation  point 

— A  y  /  / 

r  •  r  =  X  dsinecos<(>  +  y  dsinesino  +  2/3  cose 


jl(r')  =  ^lA(x^)lB(y‘')IC(2') 

^  =  “^  the  transmitter  polarization 

Putting  this  all  together  yields 


(2-5) 


Aff(r^)  = 

"^±(r) /"1a(x' )e"J^xX  dxYlB(y^  dy^^C(z^  )e"J^i®  dz*"  (2-6) 

fA(^x)  FB(^y)  EcOz) 

where;  Fa(0x5  =  fourier  transform  of  1a (x^) 

ZB(<3y)  -  fourier  transform  of  lB(y^') 

£c(^z)  =  fourier  transform  of  lc(z^) 


We  can  now  further  simplify  Eq  (2-6)  for  Aff(r  ) 


Aff(r^)  =  -^4'(r)EA(/3x)EB(^y)Ec(^z) 


(2-7) 


ft* 

*  ^  . 
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As  an  example,  Isolating  just  row  0:  m  — >  14  to  ~14,  n  =  0, 
each  isotropic  source  can  be  represented  by  a  delta  function 


(see  Figure 

10) 

-14 

row  0 

<— > 

EiaCx  ) 

m=14 

i(x^  + 

ml)] 

EB(dy) 
row  0 

<— > 

i(y'  + 

nl)| 

Fc(3z) 
row  0 

<— > 

IC(z') 1 

[4(2'^  + 

q)] 

where :  1 , q 

are  defined  on  Figure  10 

Ia  =  ImnqeJ^““> 
IB  = 

IC  =  ImnqeJ^‘l“) 


Imnq  =  amplitude  of  each  isotropic  point  source 
according  to  the  normalized  aperture 
distribution 

a  -  interelement  progressive  phase  difference  («=  0) 


Combining  all  terms  for  row  0  yields 

Aff(r')  =  -^i(r)Inmq|  L  (2-8) 


Finally,  for  row  0,  dropping  the  subscript  q  on  Imnq  since  it 
is  constant  and  n  =  0 

Iff(r')  =  -‘^i(r)Inin[  (2-9) 

m=14 


But  since  there  are  29  rows,  n  — >  14  to  -14 

Aff(r'')  = 

-^i(r)Imn|  e“J^y(nl)][  e-J^zq| 


(2-10) 


To  get  lD(e,^=0)  from  iff(r'')  for  an  H-plane  cut  (  <<*  =  0)  use 

1dW.'<*=0)  =  (■?■')  -  (2-11) 

The  radial  component  of  Aff(r^)  is  disregarded 


ED(fl.0=O)  » 

2irr 

where:  Eq  is  the  initial  E-field  amplitude 
received  by  the  dish. 

see  Appendix  A  for  actual  values  for  E©,  Q,  1,  X. 


.-jdzq 


-14  -14  .  , 

E  E 

m=14  n=14 


(2-12) 


The  next  step  is  to  derive  a  normalized  aperture  distribution 
from  which  Inm  can  be  assigned. 


Normalized  Aperture  Distribution 

The  feed  antenna  is  not  isotropic.  The  effect  of  the 
normalized  feed  radiation  pattern,  Ff(r  ),  is  related  to  the 
normalized  aperture  distribution  by  this  equation  (Stutzman  and 
Thiele,  1981:420): 


Ea(r')  =  Fi^r'  =  2f  tanl"!!  ♦  ' 

where:  r^,  4^  are  spherical  coordinates  of 

the  aperture  plane  (see  Figure  10) 


(2-13; 


f  is  the  focal  distance  of  the 
reflector  (see  Appendix  A). 


The  rear  fed  half-wave  dipole  using  a  parasitic  circular 
reflector  disc  does  not  produce  a  circularly  symmetric  pattern. 
But  the  average  of  the  E-  and  H-plane  pattern  values  can  be  used 
in  Eq  (2-13)  according  to  (Stutsman  and  Thiele,  1961:428).  This 
is  shown  by  the  dotted  line  in  Figure  11,  the  normalized  feed 
radiation  pattern  at  3.12  Ghz  for  the  rear  fed  half-wave  dipole 
using  a  parasitic  circular  reflector. 


E  — plane  pattern  value 
H— plane  pattern  value 
ave  pattern  value 


Figure  11.  Normalized  Feed  Radiation  Pattern 


To  obtain  Eg^Cr  ),  the  normalized  aperture  distribution,  r  was 
calculated  in  increments  of  =  1®  using  Eq  (2-13).  Ff(r^) 
was  read  off  Figure  11,  then  Eq  (2-13)  used  again  to  produce  a 
plot  of  Ea(r^)  (see  Figure  12  next  page).  Finally  Figure  12 
was  used  in  conjunction  with  a  scale  drawing  of  the  isotropic 
radiators  in  the  aperture  plane  to  assign  ®  function  of 

r^.  The  next  step  is  to  derive  Esi2(^  »'<*=0) . 


phygAggj  Qptiga:  Discs  Only 

The  scattering  by  a  metal  disc  will  be  derived  first. 

For  a  perfectly  conducting  body  the  physical  optics  approximate 
surface  currents  are  (Stutzman  and  Thiele,  1981:455): 


2(^  X  Hi)  in  illuminated  region 
0  in  shadow  region 


(2-14) 


Figure  13  will  be  used  as  the  coordinate  system  in  the 
development  of  Esi2(^>^=0)  1fb(^»^=0)- 


Figure  13.  Coordinate  System  For  Discs  and 
and  Feed  Blockage  Terms 
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Once  £si2(®>^=0)  and  Efb(^»^=0)  developed,  these  terns 
will  then  be  phase  referenced  back  to  the  system  coordinate 
reference,  z=0  on  Figure  10.  Assume  for  now  the  flat  metal  disc 
lies  In  the  xy  plane,  having  dimensions 

p  =  disc  radius 
p  =  0«>  to  360® 
e  -  90°  to  -90® 


The  disc  Is  Illuminated  by  a  plane  wave  from  the  transmitter 
having  fields 

Il(r)  =  ‘  (2-15) 

where:  =  angle  of  Incidence  wrt  z  axis  and  ^ 

r  =  vector  from  origin  to  observation  point 
Eqs  =  Initial  E-fleld  amplitude  received  by  disc 

Assuming  free  space 

Il(r)  =  ^  X  (r)  (2-16) 

-JwPo 

Then 

Bl(r)  =  Eqs^-^cosOi  -  -Ssln^i^e-J^CxslnJl  -  zcos^l)  (2-17) 

where:  "no  =  fa>Po 
0 

Thus  the  surface  currents  on  the  disc  are  given  by  the  physical 
optics  approximation 


since  ^  ^  and  ^  =  0  at  the  surface 


jl(r‘')  =  -y  2E2aCos9ie”J^^®^*^®^  (2-18) 

Note:  The  physical  optics  approximation  is  not  exact  and  Is 
useful  when 

P  >  X 

actual  p  =  5  inches 
actual  X  =  3.79  Inches 


Neglecting  effects  of  the  edges  on  the  surface  currents,  the 
magnetic  vector  potential  for  the  scattered  field  can  then  be 
calculated  by  (Johnson,  1984) 


Iff(r^)  =  _yy'MrtiI(r')e-J<*R  da' 

SCurce-"^-  ~  - 


(2-19) 


where :  R=lr-r|®r-^»r 

r  =  |r|  =  distance  from  origin  to  observation  point 
da'  =  p'dp'  d^^ 


Collecting  all  terms  together 

(2-20) 

•  ^  p'dp' 

where:  x  =  p  cos^ 

^  =  p^  sino  |cos0cos0^  +  sin^^sin^j 

Therefore 


Aff(r')  = 


cos5ie~J^^  sin^ig-Jdr 


Aff(r  )  =  -^2poEftsC03 9^ e~ 

^o^Jrr 

P  2T 


j J e-J^p'sinfficos^'  -  0  '  )  p'  dp*  d<>* 


(2-21) 


Because  of  the  symmetry  of  the  disc,  the  scattered  fields  will 
have  no  0  dependence,  set  ^  =  0.  Therefore,  the  two  exponential 
exponents  In  Eq  (2-21)  can  be  combined  and  simplified  to 

eJ cos  (sin  -  sin01)| 


Thus  rewritting  Eq  (2-21)  and  regrouping  it  more  conviently 

(2-22) 

T)o^ 

2r 


Aff(r')  =  -^MoEosCOS0ie~'i<**‘ 
- rj^T - 


Pr  27r  /  / 

y  eJ^P  ( 
o  L  o 


ain0  -  Sinai)  dp' 


0  L  0 

From  (Stutzman  and  Thiele,  1981:566) 


7ir 


Jo(x)  =  _x  Je^xcosa  da 


2t 


where:  Jo(x)  is  the  bessel  function  of  the  zero  order 


Letting  X  =  ^p'(sina  -  ai)  and 


Affix')  -  -^poEosCOs^ie"J^*’/'/Jo 

TJor  J 


dp (sin a  -  sinai) 


dp'  (2-23) 


From  (Stutzman  and  Thiele,  1981:566) 


xJi(x)  =  y* xJo(x)dx 

where:  Ji(x)  Is  the  bessel  function  of  the  first  order. 

Letting  X  =  ^p^(slna  -  slna^)  then  dx  =  /3(sina  -  3inajL)dp' 
therefore 


Aff  (r') 


-^MoEftsCOsaie~«><^J^  r 

»?or^(sIna  -  xJ^  (x)dx 


(2-24) 


Finally 


Aff(r^)  =  -i^2AMoEoaCostfie~J^^  Jl 


Vo^Trr 

where:  A  =  area  of  disc 


J[£(^lj[ij9j;_slri£0 

^P(sTnT"^^^"sIn?^ 


(J 


Again 


Esi( 0,^=0)  a  -JwAff(r^)  disregarding  the 

radial  component 

where:  £si(9>4>=0)  Is  the  complex  E-fleld 
scattering  off  one  disc. 


Simplifying  and  combining  all  terms  yields 

lsi(0.^=0)  «  ^.1E^ode~^g»^co3»<  /2A 

2Tr  I 

where:  x  =  0p(sln$  -  slnfl^) 


How  Is  this  modified  for  two  movable  discs  on  the  reflector 
for  ♦  =  0  l.e.  (the  H-plane  cut),  see  Figure  14  next  page? 


Figure  14.  H-plane  Cut  of  Antenna  System 

Since  the  discs  are  movable,  the  phase  or  path  length  difference 

between  li( disci)  will  different  from  Ei(disc2)  when 

ti ^  t2-  Expanding  the  picture  on  disc  1  to  develop  the 
path  length  difference  (phase)  associated  with  changing  ti, 
will  lead  to  the  phase  term  to  add  onto  Eq  (2-26)  to  account 
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r  y 


-  -  -  .  *  V'V.'V  V  •  ..■V  w". 


T  n.~{  ’»."V 


phase  difference  can  be  expressed  as  (Rudduck,  1985a) 


j/3tl  I  ^cosjii  +  cos(tfii  ^  ^)j  ^ 


But  since  there  are  two  discs  the  phase  term  must  reflect  this 


Jdtl|((cos»  +  cos(«ij^  ±  e  ))] 

j  21'  2  2  -  1 


where:  subscript  1  represents  disc  1  and 

subscript  2  represents  disc  2 


After  further  algebraic  manipulation  the  phase  term  associated 
with  Eq  (2-26)  for  both  discs  in  the  H-plane  becomes 

r  V 1  f 

2J  sin  £ii^costfii  +  cosC^il  +  ®  )j  J^ljcos^il  cosCen  +  fl  )| 

+  sin  j6t2^co3ei2  +  cos(fli2  "  ®  ))  Jdt?  [cos^j^  ■*■  cos(«i2  “  ®  ^1 
eJ/3wcos9 

where:  is  added  to  reference  the  phase  term  back 

to  the  system  coordinate  phase  reference  (see  Figure  10) 


(2-27  ) 


see  Appendix  A  for  actual  values  for 


A  good  approximation  to  feed  blockage  effects  is  to  subtract 
the  radiation  pattern  of  the  parasitic  reflector  from  the 
radiation  pattern  of  the  parabolic  reflector.  Depending  on  the 
angle,  the  feed  will  block  a  portion  of  the  radiation  causing  a 
"shadow"  on  the  reflector  surface.  The  feed  can  be  modeled  as  a 
disc  with  an  area  equal  to  that  of  the  circular  parasitic 
reflector  (see  Figure  3).  Then  for  the  E-field  scattered  off 
one  disc,  Eq  (2-26)  can  be  used  with  equal  to  0. 

=  ^JEr.-fhBe-J/3rej/3qcosg(i3)/2Ji  (gr^slng)^  (2-28) 
2irr  I  (Jrosin^ j 

where:  Tq  =  radius  of  parasitic  reflector 

B  =  area  of  parasitic  reflector 

®ofb  ■  initial  E-field  amplitude 

received  by  parasitic  reflector 

gjdqcos^r  is  added  to  reference  the  phase  back  to 
the  system  coordinate  phase  reference 
(see  Figure  10) 


see  Appendix  A  for  actual  values  for 
Eofb*  ^o*  ^ 


III.  Analysis  and  Results 


gQPPut^r  fragrama 

Three  Fortran-77  programs  (see  Appendix  C)  are  used  to 
produce  relative  power  patterns  (Program  PHASENUL),  and  field 
Intensity  amplitudes  and  phase  plots  for  the  discs  alone 
(Program  PHASCANL) ,  and  the  dish  alone  (Program  PHASEQ) . 

Prograun  PHASENUL  calculates  ^  which  is 

proportional  to  the  relative  power  pattern  of  the  antenna  system 
(see  Figures  16-27). 

Iitot(9 ^  =  IEd(®»'<'=0)  +  Es12 («.<<>  =0)  -  Efb(^.'<‘=0) I ^ 

where:  Ed(®.^=0)  =  Eq  (2-12) 

Isi2(».^=0)  =  Eq  (2-27) 

Efb(«.^=0)  =  Eq  (2-28) 

Program  PHASCANL  calculates  the  amplitude  and  phase  of  the 

field  intensity  of  the  discs  alone  which  is  proportional  to 

II.  (see  Figures  30,31). 

(U) cancellation 

lEc(»,^=0)|  =  |Isi2(».^=0)|  (3-1) 

Program  PHASEQ  calculates  the  amplitude  and  phase  of  the 

field  Intensity  of  the  dish  alone  which  is  proportional  to 

Ie  (^,<^=0)1  (see  Figures  28,29). 

(Q) quiescent 


lEQ(tf,<'=0)|  =  Ied(^.^=o)  “  Efb(®.^=o)1 


(3-2) 


'y.' 


When  l£(;(0,^=O)|  is  multiplied  by  its  phase  and 
|¥q(®.<>=0)|  is  multiplied  by  its  phase,  and  the  two  terms 
added  together,  |E«roT(®t<>=0)l  results.  IEtoT(®  ^  is 

obtained  by  squaring  the  real  part,  adding  the  square  of  the 
imaginary  part,  normalizing,  and  plotting  in  dB.  By  breaking 
|Exot(®**<*=0)|  2  into  IEc(® and  I1q(0,0=O)|  the 
phase  cancelling  in  the  side  lobes  of  the  quiescent  H-plane 
pattern  of  the  dish  by  the  discs  is  more  apparent.  This  will  be 
discussed  in  the  conclusion  section. 

The  results  of  the  computer  progreuns  and  actual  measured 
power  patterns  follow.  For  comparison  purposes,  the  measured 
power  pattern  is  shown  first  followed  by  the  theoretical  pattern 
for  the  same  disc  positions.  The  system  quiescent  power  pattern 
is  generated  with  the  two  discs  flush  on  the  reflector  surface. 
Appendix  D  shows  the  actual  disc  positions  that  were  used  to 
produce  a  measured  pattern  at  a  specific  null  angle  (e). 

The  FHASENUL  program  uses  the  actual  disc  positions  to  produce 
the  theoretical  patterns.  Following  the  power  patterns  are  the 
theoretical  field  intensity  plots.  There  are  no  actual  measured 
plots  for  comparison.  The  amplitudes  of  the  field  intensity 
plots  are  not  normalized,  to  show  that  the  cancellation  effects 
of  the  discs  occur  at  the  sldelobe  level  of  the  reflector. 
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Figure  20.  Measured  H-plane  Power  Pattern,  System  Nulls 
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Figure  28.  Theoretical  Amplitude  of  H-plane  (Dish  and  Feed 
Blockage  Only):  Quiescent  Pattern 
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Figure  29.  Theoretical  Phase  of  H-plane  (Dish  and  Feed 
Blockage  Only):  Quiescent  Pattern 


Note:  When  real  part  of  field  intensity  negative  or  when 

real  and  imaginary  part  of  field  intensity  negative 
add  180°  to  phase  angle  due  to  properties  of 
arctan  function  used  in  graph  above. 
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Figure  30,  Theoretical  Amplitude  of  H-plane  (Discs  Only) 
Cancellation  Pattern 


FAR  FIELD  ANGLE  (THETA)  (Deg) 


Figure  31.  Theoretical  Phase  of  H-plane  (Discs  Only) 
Concellation  Pattern 


Note:  When  real  part  of  field  intensity  negative  or  when 
real  and  imaginary  part  of  field  intensity  negative 
add  180®  to  phase  angle  due  to  properties  of 
arctan  function  used  in  graph  above. 


in  Eflttfirna 


The  system  quiescent  power  pattern  Figure  17  compares 
favorably  with  the  actual  pattern  Figure  16.  The  theoretical 
power  pattern  side  lobe  levels  are  down  considerably  from  the 
actual  ones.  This  Is  a  result  of  an  Idealistic  model^  l.e.  no 
Interference  sources  such  as  edge  Illumination  or  diffraction 
are  Included.  For  all  theoretical  power  patterns,  the  3dB  beam 
widths,  null  placements,  average  far-out  side  lobes,  and 
spillover  radiation  regions  agree  quite  well  with  actual. 
However,  the  first  side  lobes  are  smeared  Into  the  main  beam. 
This  occurs  because  the  phase  distribution  across  the  aperture 
Is  not  constant  and  the  model  Is  not  detailed  enough  for 
diffraction  and  double  reflection  effects  (Skolnlk,  1962:261). 
The  goal  of  the  model  Is  successful;  It  duplicates  system  nulls 
(within  +  5°  of  actual),  when  the  actual  measured  disc 
position  at  a  specific  system  null  angle  is  used  In  the  computer 
program  PHASENUL.  Cancellation  up  to  25dB,  along  with  minimum 
main  beam  degradation  and  mlnumum  distortion  of  the  overall  far 
field  pattern.  Is  evident  In  the  theoretical  patterns  matching 
the  actual  patterns.  Especially  successful  Is  Figure  21.  This 
theoretical  pattern  shows  the  symmetric  nulling  effects  of  the 
actual  pattern  Figure  20  and  predicts  the  system  nulls  exactly! 
Acccordlng  to  Jacavanco,  this  was  the  first  time  he  had  seen 
symmetric  nulls  for  a  movable  disc  system  (Jacavanco,  1985c). 

For  all  theoretical  patterns  the  intensity  or  relative  depth  of 


a  null  In  dB  (uimodlfled  minus  modified),  cannot  readily  be 
compared  to  actual.  Again  this  is  a  result  of  an  idealistic 
side  lobe  level  in  the  theoretical  patterns. 

A  requirement  of  any  nulling  technique  is  that  no  main  beam 
corruption  occur  such  as  in  Figure  26.  From  actual 
measurements,  main  beam  degradation  does  not  occur  for  this 
system  until  the  system  null  angle  exceeds  approximately  50*^. 
Looking  at  Figure  26  this  is  evident.  The  system  is  set  up  to 
produce  a  system  null  at  -60°  but  disc  placement  is  such  that 
they  overshadow  the  feed  arm  significantly  and  edge  illumination 
effects  completely  destroy  the  integrity  of  the  pattern.  The 
matching  theoretical  pattern  Figure  27,  while  not  showing 
pattern  distortion,  produces  no  nulls  and  a  very  unreasonable 
overall  pattern  for  a  parabolic  dish.  In  larger  movable  disc 
systems  that  use  a  larger  ratio  of  focal  length  to  diameter 
(aspect  ratio)  and  a  horn  feed,  main  beam  distortion  effects  are 
less  prevalent  for  system  null  angles  beyond  50°. 

As  mentioned  before,  cancellation  effects  are  more  apparent 
by  using  the  field  intensity  graphs.  From  Eq  (2-2),  the  example 
on  the  next  page  shows  how  to  compute  a  power  point  from  the 
field  Intensity  graphs,  for  a  system  null  angle  of  (0)  =  -20° 
from  Figure  21 . 


From  side 
lobe  level 
of  dish 


0 


at  =  -20° 


_  _  j*2 

q  (2-2):  yPRjQe  +  y/PRJD^ 

(2)  \/Prjq  ~  llQ(tf,^=0)|  — >  see 
yPRJQ  =  9  X  10-6 

(3)  y/PRJD  ~  |Ic(«.^=9)|  — >  see  Figure  30 

From  discs  _  at  ®  =  -20® 

\/PrJD  “  2  X  10-5 

(4)  =  phase  of  Prjq  — >  see  Figure  29 

a  218® 

(5)  ^2  -  phase  of  Prjd  — >  see  Figure  31 
<02  s  92® 

(6)  Substituting  into  Eq  (2-2): 

(9  X  io-6)eJ218°  +  (2  x  10-5)e-J62®  =  7.9  x  10-6  -  ji.e  x  10-5 


(7)  a2  +  b2  ~  power,  then  normalizing  using  the  maximum 
value  at  =  0  yields: 

5.6  X  10-5  watts 


(8)  Converting  this  to  dB  yields:  -42.5dB 

Looking  at  Figure  20  for  the  actual  system  null  at  #  =  -20®, 
-42.5dB  compares  very  well  with  a  measured  value  of  -40dB. 
Ideally  the  discs  can  be  positioned  so  >/Prjq  =v^rjd 
is  180®  out  of  phase  with  02*  However^  in  practice  this 
would  represent  "perfect"  cancellation  which  is  impossible  due 
to  secondary  effects.  In  addition,  null  depths  beyond  40dB 
descended  into  the  noise  level  of  our  pattern  recorder. 


find  £fifid  Problems 


Problems  caused  by  a  dipole  feed  are  well  known  and  are 
solved  by  using  a  horn  feed.  For  this  system  the  parasitic 
reflector  helps  to  phase  the  rearward  energy  from  the  dipole  and 
direct  it  onto  the  reflector.  The  resulting  aperture 
distribution  was  not  uniform  without  the  movable  discs  on  the 
reflector  inner  surface.  The  degree  to  which  the  non-uniform 
aperture  distribution  affected  the  phase-nulling  of  the  discs 
was  not  investigated.  The  focal  distance  was  varied,  but  had 
negligible  effects  on  cancellation.  Only  the  relative  disc 
positions  changed  for  a  relative  system  null  angle,  for  varied 
focal  lengths. 

A  second  shortcoming  of  the  dipole  feed  is  that  some  of  the 
reflected  energy  is  perpendicular  to  that  of  the  primary 
radiation.  The  cross-polarized  radiation  causes  the  anteima 
gain  to  be  reduced  and  side  lobes  with  polarization  orthogonal 
to  the  primary  polarization  are  generated.  This  energy  is 
wasted  because  the  antenna  might  not  be  designed  to  respond  to  a 
different  polarization. 

Both  effects  can  be  minimized  by  a  shallow  reflector,  one 
with  a  large  ratio  of  focal  length  to  diameter.  The  maximum 
theoretical  aperture  efficiency  obtainable  from  a 
dipole-parasitic  reflector  feed  is  69%  (Skolnik,  1970:10-6). 
Theoretical  aperture  efficiencies  obtainable  for  horn  fed 
paraboloids  are  approximately  60%  (Skolnik,  1962:273). 


Effects  Not  Modeled 

The  repeatability  of  the  measured  patterns  for  the  seone  disc 
positions  proved  to  be  a  problem.  It  is  not  uncommon  to  produce 
a  pattern  at  a  specific  null  angle,  then  find  the  final  disc 
positions  to  be  off  by  as  much  as  0.5  cm  from  one  pattern  to  the 
next.  As  another  way  to  look  at  this,  if  one  took  the  disc 
positions  recorded  in  appendix  D,  set  the  discs  according  to  the 
data,  raui  the  patterns,  the  system  nulls  would  come  out  slightly 
different  each  time  (approximately  ±2^).  This  is  due  to  a 
multitude  of  effects,  the  predominate  one  being  the  "quiet-zone" 
around  the  pedestal  is  not  constant.  Since  the  antenna  system 
is  mounted  on  top  of  the  EM  Technology  Division's  building  at 
Hanscom  AFB,  the  diffraction  and  reflection  effects  of  the 
nearby  equipment  are  constantly  changing,  as  new  and  different 
equipment  is  set  up.  To  get  accurate  patterns  this  system  would 
have  to  be  set  up  on  a  calibrated  antenna  range. 

The  thrust  of  this  study  was  to  Investigate  in  a  more 
simplified  manner  the  phase  nulling  caused  by  the  dices  in  the 
side  lobes  of  the  power  pattern  of  the  reflector.  To  more 
accurately  model  this  phenomenom  the  diffraction  pattern  of  the 
discs  should  be  included  in  a  geometrical  theory  of  diffraction 
approach  to  overcome  the  simplified  ray  tracing  approach  used 
here.  Other  factors  that  should  be  more  accurately  modeled 
are:  aperture  blockage,  errors  in  the  paraboloid  surface,  and 

finally  a  horn  fed  antenna  for  uniform  aperture  illumination. 


One  of  the  basic  problems  with  a  movable  disc  system  is  the 
amount  of  time  the  computer  takes  to  move  the  discs  to  their 
final  position  to  minimize  the  total  power  of  the  received 
signal.  Once  the  discs  reach  this  final  position  a  system  null 
has  been  placed  in  the  direction  of  the  interf erring  source.  It 
was  not  uncommon  in  the  first  systems  to  see  15  to  20  minutes  go 
by,  before  disc  optimum  positons  were  established.  Since  then 
times  of  5  minutes  have  been  established  by  efficient  computer 
programming  techniques.  It  was  hoped  at  the  outset  of  the 
research  that  the  physical  optics  approach  might  help  alleviate 
this  problem.  It  became  apparent  that  the  "speed"  problem  is 
beyond  the  scope  of  the  thesis  presented  here.  What  the 
research  does  is  provide  a  solid  foundation  upon  which  to  build 
a  geometrical  theory  of  diffraction  approach.  With  an  accurate 
model  in  hand  optimum  disc  position  could  be  predicted  rapidly. 
The  research  here  did  not  address  disc  position  or  numbers  of 
discs  (more  than  two)  on  the  reflector  surface.  This  problem 
has  already  received  a  large  amount  of  research  time  at  the  EM 
Technology  Division  at  Hanscom  AFB.  Actual  disc  placement, 
size,  and  numbers  of  discs  are  more  system  specific.  Developing 
unified  guidelines  for  a  generalized  paraboloid  would  be 
difficult. 

This  study  accomplished  what  it  set  out  to  do,  investigating 
the  phase  nulling  effects  of  the  discs  on  the  side  lobe  region 
of  the  far  field  radiation  pattern.  The  next  step  would  be  to 


generate  field  intensity  2unplitude  and  phase  plots  for  a 
multitude  of  disc  positions  at  one  system  null  angle.  The  plots 
could  be  analyzed  to  find  the  optimum  disc  positions  for  equal 
amplitudes  and  opposite  phases,  then  the  power  pattern 
generated,  and  compared  to  a  measured  power  pattern  using  that 
optimum  disc  position.  If  the  comparison  was  successful,  then 
all  optimum  disc  positions  with  their  respective  system  null 
angles  could  be  cataloged  using  this  analysis.  This  would 
require  quite  a  bit  of  time  but  would  be  valuable  iii  addressing 
the  "speed"  problem  of  positioning  the  discs. 

One  might  wonder  how  this  system  would  be  used  ouside  of  the 
laboratory.  An  urban  communication  system  which  relies  on 


reflectors  exclusively  would  profit  by  using  single  fixed  nulls 
to  reduce  specular  reflections  from  large  buildings.  NASA  which 
has  plans  for  direct  broadcast  satellites  could  use  this 
technique  for  receiving  antennas  on  user’s  roofs  to  null  out 
Interference  (Jacavanco,  1984a: 2).  In  actual  field  tests  at  the 
EM  Technology  Division,  Jacavanco  has  suppressed  a  high  power 
mobile  S-band  source  while  maintaining  a  continuous  link  with  a 


desired  low  power  signal,  fixed  in  position.  Actually,  the 
ultimate  application  will  hopefully  come  in  space  to  overcome 
the  problems  of  the  Electrostatic  Membrane  Reflector  (page  1). 

While  this  model  is  only  a  first  order  approximation  it 
appears  to  accurately  represent  the  two  disc  system.  The  linear 
array  approach  to  model  the  reflector  was  used  by  Jacavanco 
(Jacavanco,  1984a),  but  the  physical  optics  scattering  off  the 
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discs  is  unique.  The  model  breaks  dovm  the  nulling  effect  into 
two  parts.  The  power  contribution  from  the  reflector  and  feed 
verses  the  power  contribution  from  the  discs.  It  is  shown  by 
proper  disc  placement  that  the  amplitude  of  scattering  off  the 
discs  can  be  made  approximately  equal  to  the  eunplitude  in  the 
side  lobe  region  of  the  reflector  and  their  phases  close  to 
180°  out  of  phase.  A  null  in  the  far  field  power  radiation 
pattern  at  a  specific  system  null  angle  results.  The 
mathematical  model  is  implemented  in  Fortran  77-code.  The 
theoretical  patterns  are  compared  to  experimental  patterns  with 
a  high  degree  of  success. 


find  Yfiifififi 


1. 

f  - 

3.12  Ghz: 

antenna  system  frequency 

2. 

A  = 

9.62  cm: 

free  space  wavelength 

3. 

1  = 

A/4 

4. 

8  = 

2Jr/x 

5. 

r  = 

5334  cm: 

175  ft  -->  distance  to  transmitter  from 
system  phase  reference. 

(the  observation  point) 

6. 

q  = 

49.4  cm: 

see  Appendix  B 

7. 

w  = 

36.9  cm: 

see  Appendix  B 

8. 

»i  = 

26°;  see 

Appendix  B 

9. 

^o  = 

2.54  cm: 

radius  of  parasitic  reflector 

10. 

B  = 

20 . 3  cm2 ; 

area  of  parasitic  reflector 

11. 

II 

o 

Q. 

6.35  cm: 

radius  of  movable  disc 

12. 

A  = 

126.7  cm2 

:  area  of  movable  disc 

13. 

f  = 

17.3  cm: 

antenna  focal  length  (measured  from  dish 
vertex  to  parasitic  reflector) 

14. 

ti  = 

distance 
(flush  = 

of  disc  1  from  reflector  face 
.3175  cm) 

15. 

t2  = 

distance 
(flush  = 

of  disc  2  from  reflector  face 
.3175  cm) 

16. 

Radius  of  Dish 

=  32  cm 

17. 

Aperture  Plane 

Area  =  3217  cm2 

18. 

Eo-- 

normalized  to  1.0 

19. 

Eos  = 

normalized  to  0.0394  (figured  as  the  ratio  of  the 
area  of  a  disc  to  the  area  of  dish  aperture) 

20. 

Eofb' 

normalized  to  0.0063  (figured  as  the  ratio  of  the 

area  of  parasitic  reflector  to  dish  aperture  area 
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Fortr, 


M# 


Program  PHASENUL  (page  62)  calculates  in  one  degree  increments 
1Etot(^»<^=0)|  ^  for  -90°  ±6  ^  +90°.  All  programs  were 
compiled  and  executed  using  IBM  PC-DOS  2.1  on  a  Corona  Home 
Computer.  PHASENUL  will  execute  upon  typing  the  file  name 
PHASENUL  at  the  system  prompt.  The  data  file  CURRDIST.DAT  must 
reside  on  the  same  floppy  disc  as  PHASENUL.EXE.  CURRDIST.DAT 
contains  the  array  of  normalized  amplitudes  of  the  isotropic 
radiators  which  sample  the  aperture  plane  at  quarter 
wavelengths.  The  program  asks  for  E©,  Eqs,  E^fb,  ti , 
t2,  (see  Appendix  A)  on  screen  as  input  from  the  keyboard  in 
single  precision  accurracy.  PHASENUL  outputs  a  power  point  per 
degree  fl  in  a  file  called  EPV?RPT.DAT.  The  system  quiescent 
power  pattern  is  generated  by  assigning: 

Eo  =  1.0,  Eos  =  0.0394,  Eofb  =  0.0063,  ti  =  0.3175, 
t2  =  0.3175.  All  other  power  patterns  are  generated  by 
leaving  Eq,  Eqs,  Eofb  fixed  at  1.0,  0.0394,  0.0063 
respectively  and  assigning  t^  and  t2  from  Appendix  D.  The 
bessel  function  in  all  programs  is  approximated  by  its  series 
expansion  to  the  third  term. 

Program  PHASCANL  (page  68)  calculates  in  one  degree  increments 
llc(e.^=0)l  for  -90°  <  0  <  +90°. 

PHASCANL  will  execute  upon  typing  the  filename  PHASCANL  at  the 
system  prompt.  The  program  asks  for  Eos»  f-l*  ^2  (s®e 


Appendix  A)  on  screen  as  input  from  the  keyboard  in  single 
precision  accurracy.  PHASCANL  outputs  a  phase  point  and 
amplitude  point  per  degree  fl  in  files  called  PHASEC.DAT  and 
AMPL.DAT  respectively.  The  plots  were  generated  by  leaving 
Eos  “  0.0394  and  assigning  t^  and  t2  from  Appendix  D. 

Program  PHASEQ  ( page  71)  calculates  in  one  degree 
increments  |1q(0,0=O)|  for  -90°  <  0  <  +90°.  PHASEQ  will 
execute  upon  typing  the  filenaune  PHASEQ  at  the  system  prompt. 
The  program  asks  for  Eqi  Eog,  E^fb*  "tl*  ^2  (s®® 

Appendix  A)  on  screen  as  input  from  the  keyboard  in  single 
precision  accuracy.  PHASEQ  outputs  a  phase  point  and  amplitude 
point  per  degree  e  in  files  called  PHASEQ.DAT  and  AMPL.DAT 
respectively.  The  plots  were  generated  by  leaving  E©  =  1.0, 

Eos  “  Eofb  =0.0063,  tj  =  0.3175,  and  t2  =  0.3175. 

This  generates  the  quiescent  plots  without  discs. 

All  graphs  were  generated  by  using  Lotus  123  (Lotus 
Development  Corporation,  Cambridge,  Massachusetts) 


PROGRAM  PHOSENUL 


Input  variables  (real): 
Qi 
L? 

Ks 

LOMt  OR  I 
F; 

BETfiX I 
BETRYl 
BETflZ: 

THETRs 

R: 

E0: 

IXY: 

Output  variables! 

EXYi 


INTEGER  R,  B,  C,  D,  E,  F,  DEG  (  Iflt  )  ,  THETRDER,  M 

REAL  O, L, K, LAMBDA,  FREO, BETAX,  BETAZ, THETA,  R,  E0,  SINTH,  COBPLTHI 
REAL  IXY (29, 29)  ,  PI,  BESSEL,  XI,  COSTH,  COSMITHI ,  BESARGl ,  SI NTHI 
REAL  BESARG2, T1 ARG,  T2ARG.  SINTIARG, SINT2ARG, ESCONST, COSTHI 
REAL  THETAI , AREA, RO, T I , T2,  W, B1 ,  R0,  U, ENORM, EPWRPT 
REAL  BESFNT i , BESFNT2, BESFNCT, AMAG (181), RMAGATO 


COMPLEX  EXY (29,  29)  ,  COLTOTAL,  ROWTOTAL,  EFIELD(29)  ,  S,  SI 
COMPLEX  ENDTH( 181 ) , EBETAZ, JCMPLX, MJCMPLX, EFBTH ( 181 )  ,  ETOTAL 
COMPLEX  ES12TH ( 181 ) , EJKW, EJKR, EJKTl, EJKT2, EFBCONST, EX JKOCDS 


Open  input  and  output  files 

OPEN  (2, FILE-’ CURRD1ST.DAT*, status^’ UNKNOWN’  ) 
OPEN  (3, FILE-’ END. DAT’ , STATUS-’ UNKNOWN’ ) 

OPEN  ( 9, F I LE-’ ES 1 2TH.  DAT’ ,  STATUS-’ UNKNOWN’ ) 
OPEN  ( 10, FILE-’  EFBTH.  DAT’  ,  STATUS-’ UNKNOWN’ ) 
OPEN  (11, FILE-’  AMAG.  DAT’  , STATUS-’ UNKNOWN’  ) 


Iniv ial ize 


on  55  1-1,29 

DO  50  J-1,29 

EXY(I,J)-(0.  0,0.0) 
IXVd,  J)-0.0 

continue 

CONTINUE 


■  •AnT.'.V'.v-V.V.'v'.V.V.V. 


DO  e>e  i«i,£9 

COLTOTRL-  <0.  0,  0.  0) 
ROUTOTRL*  (0.  0,  0.  0) 
&0  CONTINUE 


*  Read  current  amplitude  distribution  factors 

• 

DO  100  I>1,  £9 

RERD  <£,*,  END=1000>  (1XY<I,J),  J-=l,  £9) 
MR1TE(*,*)  (rXY<I,  J), J=l,£9> 

100  CONTINUE 

• 

CLOSE  (£> 

« 

*  Initialize  constant  parameters 

« 

D=49. 4 
L«£.  41 
LRMBOR-9.  e£ 

PI=3. 141592& 

K=£*PI/LflMElDR 
JCMPLX=CMPLX (0. 0, 1. 0» 

MJCMPLX=CMPLX  <0. 0, -1. 0) 

R=5334.0 


*  Input  initial  E-field  amplitude  distribution,  E0 

« 

WRITE<*,*) 

WRITE!*,  ’  (R)’ >•  ENTER  VRLUE  FOR  E0’ 

WRITE!*,*) 

RERD  «,  E0 
WRITE  !*, 120)  E0 

120  FQRMRT!«0’,’E0  =  ’ , F6. 4 ) 

« 

« 

*  Compute  exponential  factors  for  all  rows  and  thetas 

* 

DO  500  I-l,  181 

THETRDES-I-Sl 

*  WRITE!*,  11£)  THETRDEG 

*11S  FORMRT!’  THETR  =  » , 14,  *  DEGREES*) 

THETR=  !I-91)*PI/ie0 
BETRX— K*SIN  !THETfl) 

BETflZ=-K*COS (THETR) 

EPET0Z=CEXP  (CMPLX  !0.  0,  -BETRZ*0)  ) 

* 

DO  300  fl=l,£9 

DO  £00  P=1,E9 
1*1=15-8 

EXY!R,  B)«=CEXP  (CMPLX  <0.  0,  -PETRX*M*L)  ) 
COLTOTRLwCOLTOTRL+IXY(R, B)*EXY!R,  B) 
200  CONTINUE 

ROWTOTRL*  ROWTOTRL+COLTOTRL 


COLTOTBL=  <0.  0,  0.  0> 

CONTINUE 

S-CEXP  <CMPLX  <0. 0, -K*R> ) 

S1=JCMPLX*E0*K*S/ (2*PI*R) 

ENDTH ( I ) -SI *EBETPZ*ROWTOTftL 
WRITE<3, *)  THETflDEG,  ENDTH<I) 

ROWTOTOL=  (0.  0,  0.  0) 

CONTINUE 

FORMOT  (»  FOR  THETfl  =  ’,14,*  E-FIELD  =  ’,£F1£.  5) 


*  DO  1£5,  1«>1,  181 

*  REfiD<3,*,END=ia5>  DEG < I > , EFIELD ( 1 ) 

*  WRITE(*,*>  DEGdl,  EFIELDCI) 

1£S  CONTINUE 

« 

*  E-FIELD  OF  DISKS  1  AND  2 

*  Input  initial  E-field  amplitude  distribution,  E0S 

« 

510  WR1TE<*,*) 

WRITE!*, » <fl)’ )’  ENTER  VALUE  FOR  E0S’ 

WRITE!*, *) 

READ  *,  E0S 
WRITE  !*,S20>  E0S 

520  FORMAT!' 0’ ,’ E0S  =  ',F12.9) 


*  Input  phase  distance,  T1 

* 

530  WRITE!*,*) 

WRITE!*,'  !A)' )'  ENTER  VALUE  FOR  Tl' 
WRITE!*, *) 

READ  *,  Tl 
WRITE  !*,  540)  Tl 

540  FORMAT !«0’,’T1  »  ',FG.4) 

* 

*  Input  phase  distance,  T2 

• 

550  WRITE!*,*) 

WRITE!*,’ !A)’)’  ENTER  VALUE  FOR  T2’ 
WRITE!*, *) 

READ  *,  T2 
WRITE  !*,560)  TE 

560  FORMAT !' 0’ ,’ T2  =  ’,F6. 4) 

* 

580  R0°6. 35 

W«36.  9 

THETA I =£6*P I / 1 80 
COSTH I =caS ! THETA I ) 

SINTHI-SIN!THETAI ) 

ARER=12&.  7 

EJKR=CEXP!CMPLX !0.  0,  -K*R)  ) 


ESCONST*  (EeiS*K)  /  (PI*R) 


DO  e«>e  I'l,  lai 

THETRDES«=I-9l 
THETR=THETRDEB*P 1/160 
COSTH=COS<THETfi) 

S I NTH=S I N  » THETR ) 

COSPLTH I =COS (THETR 1 +THETR> 

COSM I TH I -COS  <  THETR I -THETR ) 

BESRRGl =K*RO* ( S I NTH+S 1 NTH I ) 
PESRRG£-K*RO* (SINTH-SINTHI ) 

EJKW-CEXP (CMPLX  <0.  0,  K*M*COSTH) ) 

T 1 RRG- ( K*T 1 /£ ) *  <  COSTH I +COSPLTH1 ) 
T£RRG=  (K*T2/£)  ♦  (COSTHI+COSMITHI  ) 
EJKT1-CEXP<CMPLX<0.  0, TIRRG) ) 
EJKT2-CEXP (CHPLX  <0.  0,  T£RRG) ) 

SINT IRRG-SIN (TIRRG) 
SINT£RRG-SIN(T£RRG> 

IF  (THETRDE6.  E(3.  -26)  THEN 
BESFNTl  =  .5 

ELSE 

BESFNTl  =  BESSEL  (BESRRGD/BESRRGl 
END  IF 


IF  (THETRDEG  .EQ.  £6)  THEN 
PESFNT2  -  .5 

ELSE 

BESFNT2  -  BESSEL (BESflRG£) /BESRRG2 
END  IF 

* 

ES1£TH(  I )  —  (ESCONST*EJKR*EJKW*COSTHl«RREfl* ( (£*BEBFNT1 
+  *SINTlRRG*EJKTl)  +  (£*EiESFNT2 

+  *SINT2flRG*EJKT2))) 

WRITE (9,*)  THETRDEG,  ES12TH(1) 

♦  WRITE  (*,&01>  THETRDEG,  ESlBTHd) 

•G01  FORMRT<’  ’, ’THETR  «  ’,14,’  E-FIELD  OF  DISKS  1»2  =  ’,£F10 

• 

600  CONTINUE 

« 

*  E-FIELD  OF  FEED 

•  Input  Initial  E-fleld  amplitude  distribution,  E0FB 
« 

650  WRITE!*,*) 

WRITE!*,’  (R)’  )’  ENTER  VRLUE  FOR  E0FB’ 

WRITE<*, *) 

RERD  *,  E0FB 
WRITE  (*,660)  E0FEI 
FORMRT (’ 0’ , ’ E0FB  -  ’,Fe.  6) 


660 


WRITEdE,*) 

THETPDEB,  EPWRPT 

• 

URITE(«,Sei> 

THETPDES,  EPWRPT 

*9®1 

FORMPTI’  •,* 

THETP  =  *,14, »  POWER  PTS  OF  E-FIELD 

see 

CONTyslUE 

leee 

END 

*  BESSEL  FUNCTION  ROUTINE 

REPL  FUNCTION  BESSEL tXl) 

REPL  Xl.BESl 

• 

BES1= ( X 1 /S) - ( X 1 **3 ) / ( e**3*E)  +  <  X 1**5) / (£**5* IE) 
BESSEL'BESl 

« 


RETURN 


PROGRAM  PHASCANL 


INTEGER  A, B, C, D,  E, F,  DEG  < 181 > , THETADEG 

» 

REAL  K, LAMBDA, THETA, R,  EGS,  SINTH, COSPLTHI 

REAL  PI, BESSEL,  XI,  COSTH,  COSMITHI,  PESARGl,  SINTHI 

REAL  BESARG2, T1AR6, TgARG,  SINTIARG, SINT2ARG, ESCONST, COSTHI 

REAL  THETAI , AREA,  RQ,  T1 ,  T2,  W,  R®, PHASE 

REAL  BESFNT 1 , BESENTS,  BESFNCT,  AMPL (181), PHASEC (181) 

» 

« 

COMPLEX  ETOtSl 

COMPLEX  ESiaTHdei ) ,  EJKW,  EJKR,  EJKTl,  EJKT2 

«  Open  input  and  output  files 

OPEN  (9, FILE-’ ESiaTH.DAT*,  STATUS=’ UNKNOWN* ) 


E-FIELD  OF  DISKS  1  AND  2 

Input  initial  E-fleld  amplitude  dlatribut ion,  E8S 
510  WR1TE<«,*) 

WB1TE(*,’  (A)’ )’  ENTtR  VALUE  FOR  E0S’ 

WRITE!*,*) 

READ  *,  E0S 
WRITE  (*,520)  E0S 

520  FORMAT <• 0*, ’ E0S  -  •,F12.9) 

* 

*  Input  phase  distance,  T1 

• 

530  WRITE <♦,*) 

WRITE!*, • (A)* )’  ENTER  VALUE  FOR  Tl* 
WRITE!*,*) 

READ  *,  Tl 
WRITE  (*,540)  Tl 

5*0  format (’0*,’T1  -  *,F6.*) 

* 

*  Input  phase  distance,  T2 

550  WRITE !♦,♦) 

WRITE!*,  ’  !A)’)’  ENTER  VALUE  FOR  T2’ 
WRITE!*,*) 

READ  *,  T2 
WRITE  !*,560)  T2 

560  FORMAT !’ 0’ ,’ T2  -  ’,F6.*) 

s 

580  RO-e. 35 

w-3e. 9 

PI-3. 1*15926 
THETAI «26*PI/ 180 


COSTHl -cos  <  THETft I ) 

SlNTHI-SlN<THETfiI) 

PREP- 186.7 

LPnSDP-9. 68 

R-S334. e 

K-2*PI/LPMBDP 

EJKR-CEXP<CMPLX  <0. 0, -K*R) > 

ESCONST- (E0S«K) / (PI«R) 

PMPGPTO-. 000005696 

ESCONST- <E0S«K> / <PI*R> 

DO  600  I-l,  181 

THETPDEB-I-91 
THETP»THETPDEG*P 1/180 
COSTH-COS  <  THETP ) 

SINTH»S1N<THETB) 

COSPLTH 1 -COS  <  THETPl +THETP) 

COSM 1 TH I -COS < THETP I - THETP ) 
BESRRGl=K*RO* (SINTH+51NTHI ) 
BESPRQa-K*RO* (BINTH-SINTHI ) 

EJKW-CEXP (CMPLX (0.  0,  K*W*COSTH) ) 

TIPRG- (K*Tl/2) * (COSTHI+COSPLTHl > 
T2PRQ- (K*T8/2) * (COSTHI+COSMITHI ) 
EJKTl-CEXP (CMPLX (0. 0.  TlPRG) ) 
EJKT8-CEXP (CMPLX (0.  0,  T£PRG) ) 
SINT1PRG=SIN(T1PRG) 

S I NT2PRG-S I N ( T8PRG ) 

IF<THETPDEG.EQ.  -86 >  THEN 
BESFNTl  -  .5 

ELSE 

BESFNTl  «  BESSEL  (BESPRGD/BESPRGl 
END  IF 


IF  (THETPDE6  . EQ.  86)  THEN 
BESFNT8  -  .5 

ELSE 

BESFNT8  -  BESSEL (BESPRG8>/BESPRG8 
END  IF 

ES18TH(I)-(ESC0NST*EJKR*EJKW«C0STHI*RREP« ( (8*BESFNT1 
*SI NT 1 PRG*E JKT 1 ) ♦ (8«BESFNT8 
*SINT£PRB*EJKT8> ) ) 

URITE(9,«)  THETPDEG,  ESlETHd) 

WRITE(«,601)  THETPDEG,  ES18TH(I> 

FORMPTC  THETP  -  ’,14,’  E-FIELD  OF  DISKS  1 -  ’,8F10. 8) 
CONTINUE 


OPEN  (13,FILE»'  PMPL,  DPT’  , STPTUS-’ UNKNOWN’ ) 
OPEN  (14, FILE-’ PHPSEC. DPT’  , STPTUS-’  UNKNOWN’  > 


TOTAL  E-FIELDS 


I  DO  700  I>1,  181 

THET0DEG«I-S1 

ET0T0L=ES12TH<I> 

PHOSE»ATON<filMOG(ETOTAL) /REfiLlETOTOL) > 

PHRSEC  ( 1  >  •=  <PHRSE*360)  /  (a*PI  ) 

WRITE  ( 14,  •)  THETfiDEG,  PHOSECd) 

I  CONTINUE 

I  DO  800  I>1,  181 

THETODEG=I-91 

ETOTfiL=ESl£TH<l) 

RMPL  < I ) » ( REAL  <  ETOTRL ) **£+R I MAG ( ETOTAL )«♦£)** 
WRITE (13,*)  THETfiDEG,  RMPL (I) 

I  CONTINUE 

lO  END 

BESSEL  FUNCTION  ROUTINE 

REAL  FUNCTION  BESSEL (XI) 

REAL  X 1 , BES 1 

BESl«(Xl/2)-(Xl**3)/(£**3*£)+<Xl**S)/ (£**5*1£> 
BESSEL^BESl 


RETURN 


’  '1.  ’  '^’1  •■.  fk/'-  *.r.  v*nr  htt 


^3 

V 


i® 


PROGRPM  PHOSEO 

INTEGER  P,  B,  C,  D,  E,  F,  DEG  ( 181  )  ,  THETPDEG,  M 


REPL  O,  L,  K,  LPMBDP, FREQ, BETPX,  BETPZ,  THETP, R, EP, SINTH, COSPLTHI 
REPL  I  X Y  <  £9, 89 ) , P I , BESSEL, X 1 , COSTH,  COSN I TH I , BESPRG 1 ,  S I NTH I 
REPL  BESPRGE, TIPRG,  T£PRG, SINTIPRG, SI NT£PRG, ESCDNST, COST  HI 
REPL  THETPI,  PREP,  RO,  Tl,  T£,  W,  Bl,  RP,U,  PHPSEDdBl  > 

REPL  BESFNTl,  BESFNT2, BESFNCT, PMPL  <181 ) , PHPSE 


COMPLEX  EXY (£9, £9)  ,  COLTOTPL,  ROWTOTPL,  EF lELD (£9) , 5, SI 
COMPLEX  ENDTHl 181 )  ,  EPETP2,  JCMPLX, MJCMPLX, EFBTH ( 181) , ETOTRL 
COMPLEX  ES1£TH (181), EJKU, EJKR, EJKT 1 ,  E JKT£, EFPCONST, EX JtiQCO? 


«  Open  input  and  output  files 


ODEN 

OPEN 

OPEN 

OPEN 


<£, FILE=' CURRDIST. DPT*  ,  BTPTUS=’  UNKNOWN’  ) 
<3, F ILEe*  FND. DPT* , STPTUS=»  UNKNOWN*  ) 

(9, FILE=*  ES1£TH. DPT*  ,  STPTUS=*  UNKNOWN’  ) 

( IP, F1LE=’  EFBTH.  DPT’  ,  STPTUS=*  UNKNOWN* ) 


*  Initialize  at-rays 


5P 

55 


DO  55  1=1, £9 

DO  5P  J=l,£9 

EXVd,  J)  =  <P.P,P.P) 
IXVd,  J)=O.P 
CONTINUE 
CONTINUE 


8P 


DO  6P  1=1, £9 

COLTOTOL=  <P.  P,  P.  P) 
ROWTOTPL"  <P.  P,  P.  P) 
CONTINUE 


*  Read  current  amplitude  distt  ibut ion  factors 


IPP 


DO  IPP  1=1,  £9 

REPD  <£, », END=1PPP>  dXY<I,J),  J=l,  £9) 
WRITE!*,  *>  (IXVd,  J),  J=l,£9) 

CONTINUE 


CLOSE  <£) 

« 

*  Initialize  constant  parameters 


0=49.  A 
L=£. 41 
LPMPDR=g. 6£ 
PI=.T.  1*15926 
K=£*PI /LOMPDP 


.V 

■p 


71 


JCMPLX-CMPLX (0. 0, 1 . O) 
MJCMPLX^CMPLX  (0.  0,  -K  0) 
R«5334.  0 


*  Input  initial  E-field  amplitude  distribution,  E0 

* 

WRITEC*,*) 

WRITE <*,• (fi)’ )’  ENTER  VPLUE  FOR  E0’ 

WRITE!*, *) 

RERD  «,  E0 
WRITE  <*,1E0)  E0 

120  FORMRT!’ 0’ , ’ E0  =  •,F6.A) 


*  Compute  exponential  factors  for  all  rows  and  thetas 

• 

DO  S00  I«l,  181 

THETfiOEG=I-91 

*  WRITE!*, 112)  THETRDEG 

*112  FORMRT!’  THETR  =  ’,14,  ’  DEGREES’) 

THETR-  ! 1-91 )*PI/ 180 
BETRX»-K«SIN!THETfi) 

EETflZ=-K*COS !THETR) 

EEETRZ«=CEXP  (CmPLX  <0.  0,  -BETRZ*0)  ) 

« 

DO  300  R=l,29 

DO  200  8=1,29 
M=15-P 

EXY(0,  B)=CEXP!CmPLX  !0.  0,  -E<ETRX»M*L)  ) 
C0LT0TRL=C0LTOIflL+IXY!O,  B)*EXY(R,  E) 
200  CONTINUE 

ROWTOTRL=ROWTOTflL+COLTOTRL 
COLTOTRL=!0. 0,0.0) 

300  CONTINUE 

» 

S«CEXP !CMPLX (0. 0,  -K*R) ) 

S1=JCMPLX*F0*K*S/ !2»PI*R) 
ENDTH!I)=Sl*EDETRZ*R0WTOTflL 
WRITEI3,  *)  TMfRDEG,  ENDTH  (I ) 


ROWTOTRL=  !0.  0,  0.  0) 


500 

CONTINUE 

122 

FORNRT  !’  ’,’  FOR  THETR 

«  ’  ,  14,  ’  E-FIELD 

« 

« 

DO  I2S,  1=1,  181 

« 

READ! 3, ♦, END=125) 

DEG!1),EF1ELDH) 

« 

WRITE!*,*)  DEG<I), 

EFIELD!! ) 

125 

CONTINUE 

« 

*  E-FIEl.D  OF  DISKS  I  RND  2 

« 


Input  initial  E-field  amplitude  distribution,  E09 


5J0  WRITE(*,*) 

WRITE(*,’ (fi)’ )•  ENTER  VALUE  FOR  EOS’ 
WRITEC*, *) 

READ  *,  EOS 
WRITE  (*,520)  E0S 

5E0  FORMAT (’ 0’ E0S  =  ’,F12.  9) 


*  Input  phase  distance,  T1 

• 

530  WRITE(*, *) 

WRITE<*,  MA)*  )’  ENTER  VALUE  FOR  Tl* 
WRITE (*, *) 

READ  *,  T1 
WRITE  (*,540)  T1 

540  FORMAT (’ 0’ ,’ T1  =  ’ , F6. 4) 

* 

*  Input  phase  distance,  T£ 

* 

550  WRITE (*,*) 

WRITE(*, • (A) ’ ) ’  ENTER  VALUE  FOR  T£’ 
WRITE (*, *) 

READ  *,  T£ 

WRITE  (*,560)  T£ 

560  FORMAT <’ 0’ ,’ T£  =  *,66.4) 

* 

580  RQ=6. 35 

W=36. 9 

THETAI=£6*PI/ia0 
COSTH I =COS ( THETA  I ) 

.  SINTHI=S1N(THETAI) 

AREA=1£6. 7 

EJKR=CEXP(CMPLX<0. 0, -K«R) ) 
ESCONST=(E0S*K) / (P1*R) 

* 

DO  600  1=1,  181 

THETADE0=I-31 
THETA=THETADeG»PI/180 
COSTH=COS (THETA) 

SINTH=SIN (THETA) 

COSPLTHI=COS  (THETAI-* THETA) 
C0SMITHI=C0S(THETAI- THETA) 
BESORGl=K«RO*  (SINTH+SINTI’I ) 
BESARG£=K*RO* (SINTM-5INTHI ) 
EJKW=CEXP(CMPLX  (0.  0,  ^,»W»C03th)  ) 
T1ARG= (K*Tl/£) • (COSTHI*COGPLTHI ) 
T£ARG=  <K*T£/£) * (COGTHI+COSMI Th^ ) 
EJKT1=CEXP(CMPLX(0.  0, TIARG) ) 
EJKT£=CEXP(CMPLX (0. 0, TfPRG) t 
SINT1ARG=S1N(T lARG) 

SINT£ARG=5IN (TSARG) 

« 

IF (THETADEG. EO.  -£6)  THEN 
BESFNTl  =  .5 


BESFNTl 
END  IF 


BESSEL ( BESARG 1 ) / BESARG 1 


IF  (THETADEG  .  EO.  26)  THEN 
BESFNT2  »=  .5 

ELSE 

BESFNT2  =  BESSEL <BESARB2)/BESARG£ 

END  IF 

ES1£TH( n=- <ESCONST*EJKR*EJKW*COSTHI*AREA« ( (£«PESFNT1 
*SINT1ARG*EJKT1 )♦ r£*PESFNT2 
*SINT2ARG*EJKT£) ) ) 

WRITEO,  *)  THETADEG,  ESIETH(I) 

WRITE**, 601)  THETADEG,  ESIETHCI) 

FORMAT*’  THETA  =  ’,l<i,’  E-FIELD  OF  DISKS  1  =  ',£:F10.  8) 


CONTINUE 


E-FIELD  OF  feed 


Input  initial  E-field  amplitude  distribution,  E0FB 
WRITE**, *) 

WRITE**, ’  *A)’ )’  ENTER  VALUE  FOR  E0FB’ 
WRITE**,*) 

READ  *,  E0FB 
WRITE  **,660)  E0FB 
FORMAT*’ 0’ ,’ E0FB  =  ’  ,  FO.  6) 

Bl=20. 3 
R0=2.  54 


DO  700  1=1,  lai 

THETADeG=I-3t 

THETA=THETADEG*PI/Ifl0 

SINTH=SIN(THETA) 

C0STH=C0S< THETA) 

EFPC0NST=CMPLX (0.0,  <E0FP*K)/*e*Pl*R) ) 
EX  JKaCOS=CEXP  (CMPLX  *0.  0,  K«(?«COBTH)  ) 
U=K*R0*SINTH 


IF  (THETADEG  .  EO.  0)  THEN 
BESFNCT  =  .5 

ELSE 

BESFNCT  =  BESSEL <U)/U 
END  IF 


EFBTH  *  I )  =EFPCONST«CJI'.R*EX  JKC'CDS*P1  *£«  ( PESrNCT } 

WRITE (10,*)  THETADEG,  EFBTH *  I ) 

WRITE**,  701)  THETADEG,  E'‘BTH*I) 

FORMAT*’  ’,’ THETA  =  E-FIEI.D  OF  FLF.D  =  ',,'?ri0.(:,i 


CONTINUE 


OPEN  (13, FILE=’ OMPL.  DOT* , STOTU5=*  UNKNOWN*  ) 
OPEN  (14, FILE-’  PHOSEO.  DOT* ,STOTUG=*  UNKNOWN* ) 


TOTOL  E-FIELDS 

DO  750  1=1,  181 

THET0DEG=I-91 

ETOTOL=ENDTH ( I ) -ES 1 ETH ( 1 ) -EFBTH ( I ) 
PHOSE=OTON(01NDG(ETOTOL) /REOL(ETOTOL) ) 

PHOSEO ( I ) = (PHOSE»3e0) / <£*PI ) 

WRITE(14,*)  THETODEG,  PHOSEO ( 1 > 

CONTINUE 

DO  Bee  1=1,  181 

THET0DEG=I-91 

ETOTOL=ENDTH ( I ) -ES lETH ( 1 ) -EFBTH ( I > 

OMPL ( I ) = ( REOL (ETOTOL) *«£+01MPG (ETOTOL ) «*£>** 
WR1TE(13,*)  THETODEG,  OMPL ( I > 

CONTINUE 

e  END 

BESSEL  FUNCTION  ROUTINE 

REOL  FUNCTION  BESSEL ( X 1 ) 

REOL  X1,BES1 


BESSEL=BESl 


RETURN 


PATTERN 

Disc  1 

Dlsc2 

Figure  18,  ig 

System  Null  at: 

0  =  +350 

ti  =  5.4  cm 

t2  =  4.6  cm 

Figure  20.  21 

System  Null  at: 

9  =  -20°,  +20° 

ti  =  5.4  cm 

t2  =  5.2  cm 

Figure  Z3. 

System  Null  at: 

e  =  -350 

ti  =  4.9  cm 

t2  =  5.2  cm 

Figure  M. 

System  Null  at: 

e  =  -450 

ti  =  Flush 

(.3175  cm) 

t2  =  6.3  cm 

Figure  26.  27 
System  Null  at: 
e  =  -60° 


ti  =  7.6  cm 


t2  =  Full  Extension 
(8.7  cm) 
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scattering  off  a  flat  me'.al  disc  is  used  for  discs  and  feed 
blockage  effects.  The  array  theory  term  plus  the  leej  olockage 
term  yields  the  fi\ld  *ntens_iTy  quiescent  patterns,  phase  and 
amplitude,  prcporti^al  to  |  ^  ,^}  j  .  The  physical 

terra  for  thV  discs  yields  the  field  intensity 
cancellation  patte^rns,  phase  and  amplitude,  proportional  to 
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No  secondary  effects'  such  as  diffraction  or  edge 
iilualnation  are  considered. 

A  computer  code  was  written  to  implement  this 
approach  and  the  theoretical  patterns  produced. 

Actual  measured  patterns  are  comparec  to  the 
theoretical  patterns.  There  is  good  agreement 
between  theory  and  actual  measurements.  Finall> 
applications  of  this  antenna  system  are  presentee. 
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